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DETERMINATION OF A VEHICLE'S ABSOLUTE VELOCITY 

FIELD OF THE INVENTION 

The present invention relates generally to a method of, a system for and a 
computer program product for determining the velocity of a wheeled vehicle. 

BACKGROUND OF THE INVENTION 

Knowledge of a vehicle's speed is fundamental to the driver, but also for 
many modern control systems as anti-lock braking systems (ABS), dynamic 
stability systems, anti-spin systems and traction control systems. Also, recent 
approaches to driver safety information systems, as road friction indicators and 
sensor-free tire pressure monitoring systems require precise velocity 
information. 

There are no standard sensors in commercial cars that can measure 
absolute velocity. The speedometer of wheeled vehicles is based on the formula 

v = ar »om [ m /s], where <° is the wheel rotational speed and r " om is the nominal 
wheel radius. The rotational speed of a wheel is accurately measured by on- 
board sensors, as available for example in the antilock braking system (ABS). 
Wheel radii, however, may change with temperature, wear and even with the 
velocity to be measured due to centripetal forces which makes it impossible to 
exactly determine the wheel radius. Thus, in practice only a nominal 
approximate value can be used which leads to up to 10% error in the velocity 
measurement. This may be acceptable for the driver, who can tolerate such an 
error, but both the driver and the control and information systems listed above 
benefit from increased precision of the velocity measurement. 

Prior art describes different principles to measure the absolute velocity of a 
vehicle: 

In image correlation techniques, a camera provides an image of an area of 
the road surface beneath the vehicle, and by taking a new one which partly 



covers the same road surface area, image correlation provides an absolute 
velocity. While being very accurate and insensitive to wheel grip and 
surface, its drawback is its cost and sensitivity to dirt and damages. 
Navigation using GPS and digital street maps enables information of the 
driven path from which the absolute velocity can be computed. The price is 
still high, and there is no integrity guarantee for such systems when used in 
safety critical applications. 

Fix point reference based approaches, either requiring an on-board vision 
system or markers in or along the road and corresponding sensors in the car. 
Using a well calibrated extra free-rolling wheel. Though being used by road 
authorities, it is hardly a mass market solution. 

Sensor fusion of longitudinal accelerometer signal and wheel rotational 
speed, as disclosed in EP 1 274 613 A1. This approach promises accurate 
estimation of absolute wheel radius, and thus vehicle velocity, but requires a 
non-standard though cheap sensor. 
Another approach is based on measuring vehicle vibrations both at the 
front axle and at the rear axle of the car. The front axle and the rear axle feel 
road bumpiness in a time delayed manner. From the time difference between 
vibrations at the front axle and at the rear axle and the axle distance, the 
velocity of the vehicle can in principle be determined. In prior art, there are 
several implementations which are based on this approach. 

The German patent application DE 34 35 866 A1 discloses a system which 
uses the correlation of wheel suspension signals to determine a vehicle velocity. 
The suspension sensors measure up-down movement of the front and the rear 
axles to determine the velocity. The European patent application EP 1 014 092 
A2 describes a similar system based on the same type of sensors. A main 
drawback of these systems is that few vehicles are today equipped with the 
necessary sensors. 

The European patent publication EP 0 294 803 A2 discloses a system 
based on force sensors in a spring leg. These sensors are spaced in the driving 



direction, and the absolute velocity is computed by correlation analysis as 
above. 

The system disclosed in US patent publication US 5,301,130 is based on 
correlation analysis of data obtained by vertical acceleration sensors. 

Another type of system which is based on correlation analysis is disclosed 
in the German patent applications DE 27 51 01 2 A2 and DE 28 49 028 A2. The 
systems described therein use transducers for generating electrical signals 
reflecting the road surface unevenness. 

All these known systems have at least two of the following drawbacks: 
Most approaches require dedicated (costly) sensors. 
- Correlation based approaches uses Fast Fourier Transform (FFT) algorithms 
which are memory and computational intensive. 

Standard approaches to correlation analysis described in prior art using either 
time-domain or frequency domain (FFT) inevitably require that the velocity 
must be constant during each batch, which is hardly true in practice. Even 
small velocity variations imply that accuracy is lost. 
Furthermore, none of the above presented correlation techniques disclose a 
solution for the problem of varying vehicle velocities. Indicative of the vehicle 
velocity is the maximum peak of a correlation function whose position depends 
on the time delay between rear axle vibrations and front axle vibrations. Since 
the position of the correlation peak shifts with varying velocity, a smeared 
correlation peak is obtained from a data sample which was recorded with 
varying vehicle velocity. This peak shift deteriorates velocity measurements. 
But particularly in situations with varying velocity (breaking, acceleration, 
cornering, etc.) accurate vehicle velocity information is required by safety 
control systems. 

SUMMARY OF THE INVENTION 

A first aspect of the invention is directed to a method of determining the 
velocity of a vehicle having at least one pair of a front and a rear wheel which 



are spaced by a wheel spacing. For at least one pair of wheels, the method 
comprises determining front and rear wheel speed signals indicative of the time 
dependent behavior of the front and rear wheel speeds, respectively. The 
method further comprises correlating the front and rear wheel speed signals in 
order to determine a specific correlation feature indicative of the time delay 
between the front wheel and rear wheel speed signals. Finally, the method 
comprises determining the velocity of the vehicle based on said correlation 
feature and the wheel spacing. 

Another aspect of the invention is directed to a system for determining the 
velocity of a vehicle having at least one pair of a front and a rear wheel which 
are spaced by a wheel spacing. The system comprises wheel speed sensors 
arranged to provide front and rear wheel speed signals indicative of the time 
dependent behavior of the front and rear wheel speeds, respectively. The 
system further comprises a correlation unit arranged to correlate the front and 
rear wheel speed signals in order to determine a specific correlation feature 
indicative of the time delay between the front wheel and rear wheel speed 
signals; and to determine the velocity of the vehicle based on the wheel spacing 
and the correlation feature thus determined. 

Yet another aspect of the invention is directed to a computer program 
product including program code for carrying out a method, when executed on a 
processing system, of determining the velocity of a vehicle having at least one 
pair of a front and a rear wheel which are spaced by a wheel spacing. The 
program code is arranged to determine front and rear wheel speed signals 
indicative of the time dependent behavior of the front and rear wheel speeds, 
respectively. The program code is further arranged to correlate the front and 
rear wheel speed signals in order to determine a specific correlation feature 
indicative of the time delay between the front wheel and rear wheel speed 
signals. The program code is finally arranged to determine the velocity of the 
vehicle based on said correlation feature and the wheel spacing. 
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Other features are inherent in the methods and systems disclosed or will 
become apparent to those skilled in the art from the following detailed 
description of embodiments and its accompanying drawings. 

DESCRIPTION OF THE DRAWINGS 

Embodiments of the invention will now be described, by way of example, 
and with reference to the accompanying drawings, in which: 



Fig. 1 shows a four-wheeled car driving on a road with surface unevenness; 
Fig. 2 is a schematic view of an embodiment of a wheel speed sensor; 
Fig. 3 shows two typical wheel speed signals which were re-sampled to the 
time domain; 

is Fig. 4 shows the coherence function of the signals in Fig. 3 which is 
indicative of the frequency spectrum of the signals; 

Fig. 5 shows an estimated correlation function of the signals in Fig. 3. 

Fig. 6 shows the coherence function of the signals in Fig. 3 indicative of the 
correlation at each frequency; 
20 Fig. 7 shows the auto-correlation function of the signals in Fig. 3; 

Fig. 8 shows the absolute value (upper plot) and the phase value (lower plot) 
of the ratio of the Fourier transformed wheel speed signals of the 
front-left and rear-left wheels. A straight line is fitted to the curve in 
the lower plot whose slope is a measure for the time delay in the 
25 corresponding wheel speed signals; 

Fig. 9 shows, similarly to Fig. 3, two typical wheel speed signals which were 
re-sampled to the time domain, now with varying speed; 

Fig. 10 shows a schematic block diagram of the system for determining the 
velocity of a vehicle; 

30 Fig. 11 shows a block diagram of the method scheme of determining the 
velocity of a vehicle. 



DESCRIPTION OF THE PREFERRED EMBODIMENTS 



In general the presented method and system for determining the velocity of 
a vehicle is based on a correlation analysis of the time dependent behavior of 
wheel speed signals of a vehicle's front and rear wheels which varies with road 
bumpiness or unevenness. The cause for such variations may be any small or 
large road feature like asphalt texture, split, small stones, bumps, etc. These 
features induce these variations via the tire-road contact. The variations which 
are induced at a particular wheel result in a specific time dependent behavior of 
the respective wheel speed signal. The wheel speed signals are measured in a 
pair of a front and a rear wheel which are running in a lane so that they feel the 
same road features in a time delayed manner. Wheel speed sensors measure the 
front and rear wheel speed signals. The instantaneous wheel speed signals are 
influenced by the road features which were instantaneously passed by an 
individual wheel. A correlation analysis of the front and the left wheel speed 
signals shows specific correlation features, like a maximum in the correlation 
function, which are indicative of the time delay between the front wheel and 
rear wheel speed signals. From this time delay and the spacing between the 
front wheel and the rear wheel of the vehicle, it is possible to determine the 
velocity of the vehicle. Some aspects of this method are now described in more 
detail. 

The vehicle may be any wheeled vehicle, like cars, lorries, trucks, 
motorcycles, trains, etc. which have a front and a rear wheel in contact with 
ground. In the following embodiments, the two wheels which are running in a 
lane (a front wheel and a rear wheel) are referred to as belonging to a particular 
pair of wheels. In general, the front and the rear wheel are mounted on different 
axles but the axles are not required to be the first and the last axle of the 
vehicle. For example, the front-left and the rear-left wheel of a four wheeled car 
may constitute an appropriate wheel pair. In vehicles with more then two axles, 



two arbitrary axles of the entire set of axles may be chosen as the two wheels 
in a lane which are here denoted as front and rear wheel. The front and the rear 
wheel should preferably run in a lane so that during straight driving road 
features which are passed by the front wheel are subsequently passed by the 
rear wheel. Of course, slight displacements of the front and rear wheels with 
reference to each other are tolerable. Wheel speed variations at the two wheels 
of such a pair of wheels are the basis for the velocity determination as 
presented herein. More then one of the above defined wheel pairs may be 
included in a velocity analysis to enhance the performance of the system, but, in 
the following embodiments, the principles of the velocity determination method 
are presented with one pair of wheels only. 

The two axles respectively the two wheels are spaced by a distance which 
is in the following denoted as wheel spacing. The wheel spacing is referred to as 
parameter B in the figures and formulas. 

Fig. 1 shows an embodiment which illustrates the environment in which the 
method as presented above is implemented. A car 1 having a car body 4 is 
driving on a road 5 in forward direction v. Due to the forward motion of the car 
1 , the front wheels 2 and the back wheels 3 are rotating counter-clock wise. 
The current rotational speed of the front tire 2 is denoted as co and the current 
front wheel radius as r . Via the tire-road contact wheel speed variations are 
induced by road bumps 6. The distance between the front and rear axle is 
denoted as B . 

Fig. 10 shows a schematic block diagram of the system for determining the 
velocity of a vehicle. A front wheel speed sensor 1 2 and a rear wheel speed 
sensor 13 associated with a pair of wheels 11 provide a front wheel speed 
signal co front and a rear wheel speed signal co reart respectively, to a correlation 

unit 14. In the correlation unit 14, the front wheel speed signal co front and the 
rear wheel speed signal co rear are correlated to obtain the vehicle velocity v. 

Fig. 1 1 shows a block diagram of the method scheme of determining the 
velocity of a vehicle. The method scheme starts at the block which is denoted 
by 'START 1 . In a next block! 5, front and rear wheel speed signals are obtained 
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from front and rear wheel speed sensors. In a subsequent block 1 6, these front 
and rear wheel speed signals are correlated with each other. In a next block 17, 
a specific correlation feature indicative of a time delay t is identified in the 
correlation result. In a final block 18, the vehicle velocity v is obtained from the 
thus obtained time delay r and the known wheel spacing B of the pair of 
wheels. 

Wheel Speed Sensors 

The sensor which provides the wheel speed signal may be any common 
wheel speed sensor. Preferably, the wheel speed sensors of an antilocking 
system (ABS) are used in the embodiments since such ABS-sensors are already 
mounted in a majority of the vehicles today. Wheel speed sensors are well 
known to the person skilled in the art. They typically comprise a toothed wheel, 
a slotted disk or any otherwise featured disk which is rotating in synchronism 
with a vehicle wheel. A sensor is arranged to generate a signal each time a 
feature of the disk enters or leaves its range of observation. The sensor may be 
an optical sensor, a magnetic sensor (e.g. a HALL-sensor) or any other 
contemptible type of sensor. The sensor produces electrical signals which are 
transported by wires or radio transmission to a subsequent unit for further 
processing. 

Fig. 2 shows a schematic view of an embodiment of a wheel speed sensor. 
The wheel speed sensor 7 comprises a toothed wheel 8. The toothed wheel has 
L teeth 10 which are equidistantly spaced over the circumference of the wheel. 
In the embodiment as shown, there are L-l teeth which are spaced by gaps of 
the same dimensions as the teeth. Thereby, the circumference is separated in L 
identical sections which all cover the same angle cc=2/t/L. A sensor 9 is 
arranged to observe passing teeth. The sensor generates a high level signal 
when covered by a tooth and a low signal when not. Thus it generates a 
periodic pulse train having a periodicity which depends on the rotational speed 
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of the rotating disk. Preferably, a signal is generated at each flank of the pulse 
train indicating an advance of the wheel by an angle of oc=2tl/L. 

Embodiment of a Signal Model 

5 

In this section, the principles of the velocity measurement are explained 
with reference to a continuous time wheel angular speed signal oo(t). With 
Continuous time 1 signal a real or hypothetical signal is denoted which provides a 
data point for every time instance in the available time interval. Some sensors, 

10 however, provide sampled measurement signals which are not continuous time 
signals but discrete signals since, for example, in digital measurements only a 
limited number of measurement values is recorded. A discussion of an 
implementation of the velocity measurement in connection with discrete signals 
will be provided further below. 

15 Let co t {t) denote the velocity of each wheel 2 = 1,2,... of a vehicle. The 

embodiment which is shown in Fig. 1 applies to a four-wheeled car where the 
front-left wheel is numbered 1 , the front-right 2, the rear-left 3 and the rear-right 
4. Examples of wheel speed signals are shown in Fig. 3. The plot shows 
measured wheel speeds in rad/s of the four wheels of a car as function of time. 

20 The data was recorded over a time interval of 60 seconds. The wheel speeds 
are shown in the interval from 41.9 rad/s to 42.7 rad/s which illustrates that the 
fluctuations of the wheel speeds around their median are in the range of some 
percent. 40 rad/s corresponds roughly to 40 km/s which is about 1 1 m/s. Plotted 
are four curves corresponding to the four wheels FL = front-left, FR = front- 

25 right, RL = rear-left, RR = rear-right of the car. One curve contains an offset 
that is due to different tire radii, wheel slip or cornering. Both curves contain a 
barely visible disturbance that is delayed by 0.3 seconds. 

The general principle applied in the embodiments is that the wheel speed 
signals contain a road induced disturbance which appears in the rear axle 
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signals cc\(t) and <^(/) r seconds later than in the front axle signals co^it) and 
co A {t) . An embodiment for a signal model thus is 

<q(0=*(0 + rf(0 + ei(0 (D 
<%(t) = s(t) + d(t-T) + e,(t) (2) 

where s{i) is the deterministic velocity component which stems from the 
velocity variations of the vehicle during driving and typically lies in the interval 
[0-5] Hz, d(t) is a road induced disturbance caused by road bumpiness which 
typically lies in the interval [30-60] Hz and e ( (t) are all other kind of 
disturbances caused by wheel suspension, sensor errors and within the 
driveline. e.(t) may also contain small offsets caused by wheel slip on the driven 

wheels and road curvature (the front and rear wheels do not follow exactly the 
same path). This offset is clearly visible in Fig. 3. 

The spectra shown in Fig. 4 and the auto-correlation sequences shown in 
Fig. 7 are obtained after high-pass filtering the wheel speed signals to filter out 
s(t). The resemblance of the curves indicates that the common term d(t) 
dominates the noise term 

By determining the time delay t between the disturbance component d(i) 
of the front and rear wheels the absolute velocity v of the vehicle can be 
calculated from the relation r= ^/ where B is the spacing between the axles 

of the front and the rear wheel. In the embodiments as shown, this time delay t 
is determined by a correlation analysis. 

Velocity Determination via Correlation Analysis 

The cross correlation between the front and rear wheel speeds is defined 

as 
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R, 3 (t) = E[(cq (f) - Efa (0)Xoi (t - r) - - T)))] . 



(3) 



with E(cc(t)) denoting the expectation value of cc(f) ■ Small disturbances injected 
by uneven road surface will occur first on a\(t) and B/v l seconds later on 
Here v, denotes the velocity of the left side of the car. The cross correlation 
function is a function of time delay t and will show a peak at 



When driving with constant speed, we can directly obtain R u (r) • The plot of Fig. 
4 shows the correlation function R u (r) of the wheel speed signals of the front- 
left and the rear-left wheel and the correlation function R 24 (r) of the wheel 



the plot represents the time delay t in the interval from -1 to 1 seconds. The 
ordinate denotes the corresponding values of the cross correlation function 
R n (r) . Each correlation function shows a peak of maximal oscillation for a 

particular time delay value r. These peaks indicate that disturbances in the 
wheel speed signals are most similar at these particular time delay values. 
The time delay is obtained from i^ 3 (r) by 



Inserting the obtained t in equation (4) then yields the absolute velocity of the 
left side of the car and the wheel radius of the rear-left. The above presented 
method however is only applicable when the vehicle is driving with constant 
velocity during the time interval in which the wheel speed signals which are 
inputted in equation (3) are recorded. A varying vehicle velocity would result in 




(4) 



speed signals of the front-right and the rear-right wheel of a car. The abscissa of 




(5) 
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a smearing of the peak in the cross correlation function since the peak shifts 
with increasing velocity to the left and with decreasing velocity to the right. 

Correlation Analysis for Varying Vehicle Velocity 

In the following an embodiment is presented which remedies the problem 
of the shifting cross correlation peak at non-constant vehicle velocities. Since 
the correlation feature in equation (3) moves with velocity, a better approach is 
to study a velocity compensated correlation function. The idea here is to use the 
information about the vehicle velocity which is inherent in the wheel speeds. 
The wheel speed provides information about the current vehicle velocity when 
multiplied with the wheel radius, which however is an unknown scale factor. 
This information can be exploited to interchange the shifting of the correlation 
maximum with velocity by a shifting with the wheel radius. Using the relation 
x=B/cqr. and defining x^B/^ in equation (3) gives the velocity compensated 

correlation function 

7? 13 (x) = E{cq (0 - E{cq (0))(^ (t - x/cq (/)) - £(cq (t - x/cq (*)))) ■ (6) 

By maximizing ^, 3 (x)with respect to the variable x = B/r l one can directly obtain 
an estimation value r, for the wheel radius of the front-left wheel: 

g 

f } = argmax R ]3 (—). (7) 
T r 

Above, the correlation function was evaluated at r=x/cq(t) using the 
front-left wheel speed a\(t). The correlation function can analogously be 
obtained by evaluation at r=x/cq i (t) using the rear-left wheel speed a^(0 with 
x = B/r z and 

7? 31 (x) = £fo (t + x/oi (0) - £(a}(f + *M(0)fe (0 - Efa(?))) (8) 
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which gives 

g 

r 3 = argmax i? 31 (— ) - (9) 
r r 

Proceeding with the right side of the vehicle in a similar way yields f 2 and f A . 

Since it is known that the wheel radius does not deviate very much from 
the nominal wheel radius r nom it is sufficient to estimate R n (x) on a grid close 
around B/r nom . The finer the computed grid is, the better is the potential 
accuracy in determining f t . 

In summary, this embodiment allows to estimate all wheel radii r t of the 

vehicle even at rapidly varying vehicle velocities. From these wheel radii r., the 

absolute velocity v. =0%^ at each wheel can be determined and simple 

geometrical transformations can be used to find the velocity at any position of 
the vehicle. 

Resolution Enhancement 

The resolution in determining the position of the peak in the correlation 
function determines the resolution of the final wheel radius measurement 
respectively the velocity measurement. In reality, the functions cq(t) are not 
continuous time functions. Instead, they are discrete functions and their time 
resolution is mainly determined by the cog spacing of the cog-wheel in the ABS- 
sensor. Since a typical axle spacing of a car is 5 = 3111 and a typical wheel 
radius is roughly r = 0.3m the peak of R n (x) typically occurs for approximately 
x=B/r nom ~3/03=l0 which corresponds to approximately 5£/(27zr)« 3-50/2 = 75 
cogs for an ABS sensor with a cog-wheel with L = 50 cogs. In principle, the 
resolution which is obtainable when r is determined is limited by this cog 
separation. Estimating for example 76 instead of 75 cogs for the correlation 
peak position results in an radius error of approximately 1.3%. In the following, 
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embodiments are discussed which provide a finer resolution of the peak position 
determination. Another inherent problem is wheel slip on the driven wheels, 
which may introduce errors of about 1 %. For two-wheel driven wheels, this can 
be avoided by only using free-rolling wheels. 

In one embodiment, the wheel speed sensor signal is triggered both at 
positive and negative flanks of the toothed wheel. This doubles the precision 
compared to triggering only at positive flanks. 

In another embodiment which is demonstrated later, non-integer 
interpolation is used when searching for the correlation peak. 

Frequency Domain Analysis 

In the following embodiment a frequency domain approach is presented for 
determining the time delay r. It is well known to a person skilled in the art that 
a shift in the time domain corresponds to a phase change in the frequency 
domain. Fourier transforming the wheel speed signals with standard FFT- 

methods results in corresponding Fourier coefficients A time shift of t of 

the wheel speed signals in the time domain ca(i) ->co(t — T) corresponds to a 
phase shift 

&(/)-> ntfle-"** (10) 

in the frequency domain. 

One approach is to adjust a straight line to the phase curve ofQ t ^f)/Q, 3 (f) 

which should according to equation (10) have a slope 2ttt. Fig. 8 shows an 
example which illustrates this procedure. The upper plot of Fig. 8 shows the 
absolute value abs(£2,(/)/Q 3 (/)); and the lower plot shows the phase angle 
arg(Qj(/)/£2 3 (/)) of the Fourier coefficient ratio of the front-left and rear-left 

wheel speeds. A straight line is fitted to the high frequency end of the curve. 
The slope of this specific curve yields a value of r= 0.13 - 
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Another method to determine r in the Fourier domain is described in Knapp 
and Carter, "The generalized correlation method for estimation of time delay", 
IEEE Trans on Acoustics, Speech and Signal Processing, vol 24(4), pp. 320-326, 
1976, the disclosure of which is herein incorporated by reference. This method 
is based on the formula 



where W{f) is a weighting function. This weighting function can be chosen 
such that the maximum likelihood estimate of time delay is obtained. 

Time Domain and Angle Domain Signals 

In digital signal processing it is very common that signals or data are 
sampled in the time domain. In this case, measurement values like signal values 
or signal changes are stored for subsequent points in time which are 
equidistantly spaced over the time axes. However, wheel speed sensors do not 
provide their data in a way which allows a direct storage in time domain. As 
described above, a wheel speed sensor triggers a signal each time the observed 
wheel has rotated by a particular angle (oc=27i/L in the above embodiment). 
Consequently, the signal values (angles) are here equidistantly distributed over 
the signal axes, whereas the corresponding time instances are not equidistantly 
distributed over the time axes. The time intervals between two sensor signals 
depend on the rotational velocity of the observed wheel. Such type of data 
which is for example generated by sampling the output from a wheel speed 
sensor is generally referred to as event domain sampled (here angle domain 
sampled). Standard signal processing software or hardware products are 
normally not adapted to event domain sampled data. Fast Fourier Transform 
(FFT) algorithms for example are in general adapted to time domain data. 




(11) 



16 



Interpolation and Filtering 

Some of the embodiments in the previous sections are based on a 
correlation analysis in the time domain. According to the above section, the 
wheel speed signals which are originally sampled in the angle domain therefore 
have to be first transformed to time domain. A standard problem which thereby 
appears are aliasing effects. The person skilled in the art knows that aliasing 
effects occur when continuous time signals comprise frequencies above the 
Nyquist frequency of the applied sampling frequency. Such high frequency 
components of the signal may be reflected into the signal range of interest by 
aliasing effects which leads to a deterioration of the signal and the velocity 
determination performed thereon. The aliasing effects can be avoided by 
appropriately filtering and interpolating the signal when transforming from angle 
to time domain. Details about the problem of aliasing effects and respective 
remedies in the context of angle domain wheel speed signals can be found in the 
PCT patent application PCT/EP02/1 2409. Therein disclosed are methods for 
effectively transforming angle domain wheel speed signals to the time domain 
without introducing aliasing effects. The disclosure of this patent application 
shall be incorporated herein by reference since the following embodiments make 
use of the methods disclosed therein. 

The angle domain data sample is obtained by determining the time 
instances of the successive flank changes in the sensor and storing the obtained 
time instances in a measurement vector t k . Each time instance corresponds to a 

constant rotation of the wheel by an angle of =2n/L (see above section). 

The wheel speeds in the time intervals [t k _ l9 t k ] are then approximately given by 
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Manufacturing errors e. in the tooth widths of the wheel speed sensor 
(q; = 2jv/L + £ i ) lead for example to important aliasing effects which have to be 
compensated for. 

As stated above, the wheel speed signal has to be appropriately filtered. To 
a hypothetical continuous time wheel speed signal, a time continuous anti-alias 
or band-pass filter h(t) could be applied: 



In the following embodiments, equation (12) is approximated in one of the 
implementations listed below (see PCT/EP02/1 2409 for details): 

1 . Piece-wise constant co and piece-wise constant h : 




(13) 



m 



(14) 



Inserting equation (12) in equation (14) gives 



(14b) 



2. Piece-wise constant co and continuous time h : 




(15) 



This sum can be computed from tabulated values of the integral 
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3. Piece-wise linear co, with continuous time h : 

tf(0~ Zf [ ^C^_, ) + ^^-to^ ) - )) - ^Ktt (16) 

Similar to above, this can be computed from tabulated values of the integral of 
the impulse responses H(t) = £a(z)</t and H\t) = £-rt(r)rfr. 

In principle any polynomial spline function may be defined as approximation 
for co. 

The filter h may be any appropriate filter function. One embodiment is for 
example a second order Butterworth filter as used in PCT/EP02/1 2409. 

Implementation 

The correlation function of equation (6) contains expectation values which 
have to be approximated. In one embodiment of the invention, this is done by 
first high-pass filtering each wheel angular speed signal and then averaging the 
filtered output signal with an exponential filter: 

Qjr p (t) = H(z)of(t) 

In the first equation, H(z) is an appropriately designed high-pass filter, here 
expressed in the Z-transform operator z. This high-pass filter should remove low- 
frequency components in the wheel speed that come from wheel suspension 
and driver commands. If these components are included, a very strong 
correlation of all wheel signals will be found for zero lag x=0. On the other 
hand, a low-pass filtered wheel speed signal, which may be computed as 
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eqf (f) = of (O-of'CO-JiCO f approximates the signal part * (r) in (1), and is in 
(17) used for computing the normalized lag. 

In the second equation, A is a forgetting factor (the same as the pole of a 
first order filter H(z) = Xl{z - A) ) chosen to get the appropriate time constant for 
tracking a time-varying wheel radius. The forgetting factor can for example be 
chosen as A= 0.995. The purpose hereof is to recursively compute an estimate 
of the correlation function. 

Fig. 5 shows an off-line example of an estimated correlation function using 
this algorithm. There is a clear peak for a lag of 70 samples, which occurs for 
t-0.27. The ringing in the correlation comes from the auto-correlation in the 
signal itself (the signals s and J in equations (1) and (2)). In comparison thereto, 
Fig. 7 shows a similar example of an estimated auto-correlation function. As can 
be seen, the correlation function is mainly a time shifted version of the auto- 
correlation function, and much higher resolution in r can be achieved by a 
second stage where the auto-correlation is continuously (not using integer shifts) 
shifted until the best match is achieved. 

Further, the filter H(z) in (17) should attenuate all frequencies that do not 
contribute to the correlation of interest. As can be seen in the coherence 
function in Fig. 4, the frequencies around 50 Hz seem to be best correlated. 
That is, H(z) is better chosen as a band-pass filter centered around 50 Hz. 

Correlation Analysis in the Angle Domain 

Instead of transforming the wheel speed data from the angle domain to the 
time domain, the following embodiment is based on a direct computation of the 
correlation in angle domain. 

If the covariance Cov(x,y) = E(xy)- E(x)E(y) is correlated, then also 
Cov(f(x),g(y)) = E(f(x)g(y)) - E(f(x))E(g(y)) is correlated. Now, take 
/(*) = £(*) = 1/* • That is, instead of assuming tire correlation in the unit [rad/s], 
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possible correlations are investigated in the inverse unit [s/rad]. Substitution in 
equation (6) yields the modified correlation function 



R l3 (x)=E 



-J. 



1 



tq(t ) 



i 



a±(t ) 



-(^J* "'L.M'l "t* =^ -tll = 0,1,2,... 



2tt 



(18) 
-, and the 



The approximation follows from (12), where o)fc) = - 

normalized lag x is implicitly defined by *, 3 using the relation = That 

is, the integer values of / are used to update the grid points on the x-axis for the 
correlation function. 

Here, t l m denotes the time measurement at which the m-th cog of the ABS- 
sensor passes the sensor in wheel i. The function J? 13 (» is preferably 
represented by a table with a fine grid around x = B/r 3 , and the grid point 
x=B/r 3 is updated which is closest to the value t\ An embodiment for an 

algorithm in which the mean is approximated by the average cog time during the 
last complete revolution is given by 



R ]3 (x) = te n W + (1 - A)^J (fc - )- (tl +L/M -tl L/2 )) 



(19) 



Here, A again is a forgetting factor. Determining the maximum of R n (x) 
gives a good approximation of the mean inverse velocity which is insensitive to 
constant acceleration and wheel imbalance and cog errors. The grid size of the 
table containing cqr = t\ -tf should be chosen according to the requirements on 
resolution in the computed wheel radius 
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System and Computer Program Product 

The embodiments of the computer program products with program code for 
performing the described methods include any machine-readable medium that is 
capable of storing or encoding the program code. The term "machine-readable 
medium" shall accordingly be taken to include, but not to be limited to, solid 
state memories, optical and magnetic storage media, and carrier wave signals. 
The program code may be machine code or another code which can be 
converted into machine code by compilation and/or interpretation, such as 
source code in a high-level programming language, such as C++, or in any 
other suitable imperative or functional programming language, or virtual-machine 
code. The computer program product may comprise a data carrier provided with 
the program code or other means devised to control or direct a data processing 
apparatus to perform the method in accordance with the description. A data 
processing apparatus running the method typically includes a central processing 
unit, data storage means and an l/O-interface for signals or parameter values. 

Thus, a general purpose of the disclosed embodiments is to provide 
improved methods and products which enable to more accurately determine a 
vehicle's velocity by means of wheel speed sensors which are in particular 
already existing within common vehicle electronic systems (antilock braking 
system and the like). 

All publications and existing systems mentioned in this specification are 
herein incorporated by reference. 

Although certain methods and products constructed in accordance with the 
teachings of the invention have been described herein, the scope of coverage of 
this patent is not limited thereto. On the contrary, this patent covers all 
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embodiments of the teachings of the invention fairly falling within the scope of 
the appended claims either literally or under the doctrine of equivalents. 



